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ABSTRACT 



Neutral sodium is a minority species in the atmospheres of late-type stars, and line formation in local thermodynamic equilibrium 
(LTE) is often a poor assumption, in particular for strong lines. We present an extensive grid of non-LTE calculations for several 
Nal lines in cool stellar atmospheres, including metal-rich and metal-poor dwarfs and giants. For the first time, we constructed a 
Na model atom that incorporates accurate quantum mechanical calculations for collisional excitation and ionisation by electrons as 
well as collisional excitation and charge exchange reactions with neutral hydrogen. Similar to Li I, the new rates for hydrogen impact 
excitation do not affect the statistical equilibrium calculations, while charge exchange reactions have a small but non-negligible 
influence. The presented LTE and non-LTE curves-of-growth can be interpolated to obtain non-LTE abundances and abundance 
corrections for arbitrary stellar parameter combinations and line strengths. The typical corrections for weak lines are —0.1... - 0.2 dex, 
whereas saturated lines may overestimate the abundance in LTE by more than 0.5 dex. The non-LTE Na abundances appear very 
robust with respect to uncertainties in the input collisional data. 
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1. Introduction 

Sodium has established itself as an important tracer of Galactic 
chemical evolution, and numerous investigations of the Na 
abundances of late-type stars, residing in different reg i ons of 
the Galaxy have been conducted ( e.g. [T akeda et ail l2003t 
iGehren et all l2006t lAndrievskv et all I2007I) . Na is mainly syn- 
thesised during hydrostatic carbon-buming in massive stars, 
through the reaction 12 C( 12 C,p) 23 Na. As pointed out by 
IWooslev & Weaver! (Il995l) . the production is dependent on the 
available neutron excess through secondary reactions, which 
implies metal-dependent yields. In addition, there is a pro- 
duction channel via proton-captur e reactions, 22 Ne(p,y) 23 Na 
dDenisenkov & Denisenkov 1 fl990h . The latter, the so-called 
NeNa-cycle, occurs when temperatures are high enough for H- 
burning through the CNO-cycle, e.g. in the cores or H-burning 
shells of intermediate mass and massive stars. 

Abundance studies of late-type stars in the thin disk show 
a steady increase from solar to positive [Na/Fe] ratios at 
super-solar metallicities, while thin and thick disk stars instead 
form a decreasing trend slightly below solar metallicitie s (e.g . 
Edvardsson et allll993LlReddv et alll200llBensbv et all 12003: 
Shi et al., 2004). Relying exclusively on weak lines for the anal- 



ysis, LTE has been proved a reasonable approximation in this 
metallicity regime. For metal-poor stars the situation is dif- 
ferent, especially in cases where the strong Nal D resonance 
lines (588.9/589.5 nm) are the only av aila ble abundance indica- 
tors. As shown bv lTakeda etaP (120031) and lGehren et all (12006). 
[Na/Fe] ratios are slightly subsolar (—0.1... - 0.5) in metal-poor 
stars in the thick disk and the halo in the metallicity range 
[Fe/H]= -3.0... - 1.0. This deficiency is only recovered through 



non-LTE analysis, because LTE investigations tend to overesti- 
mate the abundances, sometimes by as much as ~ 0.5 dex (see 
Sect. 13. 11 1. Moreover, almost solar values are obtained from non- 
LTE analysis of extremely metal-poor stars below [Fe/H] < 
-3.0, where LTE analysis , at least of giants, i nstead yields pos- 
itive r atios, rNa/Fel ^ 0.3 (ICavrel et all 12004 lAndrievskv et all 
120071) . Interestingly. iNissen & Schusterl ( 201 Oh found evidence 
for systematic Na abundance differences of the order of 0.2 dex 
between a-poor and cr-rich halo stars, with important implica- 
tions for the presumably separate origin of these two Galactic 
components. To place all Galactic stellar populations on an ab- 
solute Na abundance scale to the same and better precision, non- 
LTE is clearly required. 

In globular clusters, Na is of particular interest, because 
the large overabundances of this element, compared with those 
of field stars of similar metallicities, imply a chemical evolu- 
tion scenario that is specific to these dense stellar systems (e.g. 
iGratton et all 12004). By detailed mapping of the Na abundance 
and its correlating behaviour with similar-mass and lighter el- 
ements we may distinguish between stars formed in different 
formation episodes in globular clusters and eventually identify 
the elusive self-enrichment process that so efficiently polluted 
the star-forming gas with the nucleosynthesis products of hot 
H-burning through the CNO-cycle and the related NeNa- and 
MgAl-chains (i.e. enhan cement of N, Na, and Al, and depletion 
of C, O, and Mg, see e.g. lLind et alJsubmittedl and lCarretta et all 
120091) . 

In this study we present ID, non-LTE calculations for several 
neutral Na lines for a large stellar grid. The given non-LTE abun- 
dances can be interpolated to arbitrary stellar-parameter combi- 
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nations and will consequently be useful for Na abundance analy- 
ses with a variety of applications. In a forthcoming paper we will 
apply the non-LTE modelling procedure described here to indi- 
vidual stars, and extend it to using 3D hydrodynamical model 
atmospheres. 

Typically, the largest uncertainties affecting the non-LTE cal- 
culations in stellar atmospheres are collisional cross -sections, 
espec ially for collisions with hydrogen atoms (e.g. lAsplundL 
2005). As described in Sect. 12.31 and 12.41 recent quantum me- 
chanical calculations of cross-sections for collisions with both 
electrons and hydrogen have significantly improved the situation 
for sodium, and a reliable atom can now be constructed. In this 
study, we put extra emphasis on assessing to which extent the re- 
maining uncertainties in collisional rates influence the statistical 
equilibrium of Na I. 



2. Non-LTE modelling procedure 

We used the code MULTI, version 2.3 dCarlssonl [19861 fl992h 
to simultaneously solve the statistical equilibrium and radiative 
transfer problems in a plane-parallel stellar atmosphere. Na was 
considered a trace element, i.e. we neglected feedback effects 
from changes in its level populations on the atmospheric struc- 
ture. The LTE assumption was used for all other species in the 
computation of background continuum and line opacity. A si- 
multaneously computed line-blanketed radiation field was there- 
fore used in the calculation of photoionisation rates. 

A grid of 764 ID, LTE, opacity- sampling MARCS model 
atmospheres dGustafsson et al.L 120081) was used in the analy- 
sis. The models span T eS = 4000... 8000 K, logg = 1.0.. .5.0, 
and [Fe/H] = -5.0... + 0.5. Sodium abundances vary from 
[Na/Fe] = -2.0... + 2.0, in steps of 0.25 dex. The highest ef- 
fective temperature is 5500 K for models with logg = 1.0, 
6500 K for logg = 2.0, and 8000 K for logg > 3.0. For mod- 
els with logg > 3.0, we adopted microturbulence parameters 
£ t = 1-0 and 2.0kms _1 and for models with \ogg < 3.0, we 
adopted £ t = 1-0,2.0, and 5.0kms _1 . All models have a stan- 
da rd composition , i.e. w ith scaled solar abundances according 
to iGrevesse et alJ ([2007), plus 0.4 dex enhancement of alpha- 
elements in all metal-poor models ([Fe/H] < -1.0). 

We define a non-LTE correction for each abundance point 
as the difference between the LTE sodium abundance and the 
non-LTE abundance that corresponds to the same equivalent 
width. Corrections are given for equivalent widths in the range 
0.01 - 100 pm. The equivalent width was obtained by numeri- 
cal integration over the line profile, considering a spectral region 
that extends +0.3 - 3 nm from the line centre, depending on the 
typical line strength. 

2.1. Structure of the model atom 

The model atom we constructed for Na consists of 22 energy lev- 
els of Na I, plus the Na II continuum (see Fig.[TJ. The energy lev- 
els are coupled by radiative transitions and by collisional transi- 
tions with electrons and neutral hydrogen. Experime ntally mea- 
sured energies were taken from the compilation by ISansonettil 
(2008) (for highly excited states isoelectronic fitting was used). 
Because we are not concerned with modelling the detailed struc- 
ture of highly excited levels, we collapsed all sublevels for 
n = 7 — 12 into super-levels. The fine-structure components of 
the 3p level were accounted for in the statistical equilibrium cal- 
culations by treating 3pi/2 and 3p3/2 as separate, collisionally 
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Fig. 1. Schematic term diagram of the 23-level Na I model atom. 
The dashed line marks the Nail continuum. The highest level 
considered in Nal is n = 12. 

coupled levelfl In addition, the fine-structure components of the 
3p-3d and 3p^4d transitions, as well as the hyper-fine structure 
components of the 3s-3pi/2 and 3s-3p3/2 transitions, were ac- 
counted for by computing the line profile as a linear combination 
of the subcomponents (see Table[T]l. 

The collisional transition probabilities were computed for 
the 3p level, not for its fine-structure components. To calculate 
these we assume that the ionisation and excitation rates of the 
sublevels were equal to that of the collapsed level. We further 
assume that the de-excitation rates to 3s scale with the statistical 
weights, i.e. the de-excitation rate from 3p3/2 is twice as large 
as from 3pi/2, and the sum of the rates is equal to that of the 
collapsed level. The collisional cross-section for electron impact 
excitation between the two fine-struct ure leve l s was estimated 
with the impact parameter method of ISeatoril (0~962). We note 
that the detailed rates are not important for the non-LTE prob- 
lem. 

2.2. Radiative transitions 

In total, 166 allowed bound-bound radiative transitions were in- 
cluded, adopting where possible oscillator strengths from the 
ab initio calculations of C. Froese Fischer Q. For most re- 
maining transitions, we used the calculated transition probabil- 
ities by K.T Taylor, as part of the Opacity Project (TOPbase 
ICunto & MendozaL 1 19921) 0. The two sets of /-values typically 
agree within 3% for the strongest, most important transitions. 
For the NaD lines, accurate experimental data exist, and we 
adopted the values listed in the NIST data bas^3 (see TableQ]). 

Photoionisation cross-sections for levels with I < 4 were 
drawn from the TOP-base (computations by K. T. Taylor). For 
the highly excited collapsed levels we adopted hydrogenic cross- 
sections. 



1 Accounting for the fine-structure components of the 3p level sim- 
plifies the numerical procedure, while als o predicting correct line for- 
mation depths for the individual lines (e.g. Mashonkina et aUl2000l) . 

2 Multi-configuration Hartree-Fock computations (MCHF). 
http : / /www . vuse . vanderbilt . edu/- cf f /mchf_collection/ 

3 http : //legacy . gsf c . nasa . gov/topbase 

4 |http: //physics .nis~ gov/PhysRefData/ASD/index .html| 
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(o) Ail is the difference in energy between the given subcomponent and the av erage energy gap brid ged by the transition. 
The energy separations of the hyperfine components of 3s-3p are taken from Svdorvk et al. (2008) 
*' T is the natural broadening width (FWHM). 

{c) cr is the broadening cross-section for collisions with neutral hydrogen at relative velocity v = 10 4 m s -1 . 
(rf) a is the velocity dependence of cr. 
w C 4 in [cmV] 



2.3. Collisions with electrons 

Cross-sections for collisional excitation by electrons can be es- 
timated using general recipes based on the Bom approxima- 
tion, which is known to overestim ate the cross-sections at low 
impact energies (e.g iParkt 1 1 97 lb - Those near-threshold ener- 
gies are most relevant for stellar atmospheres that host elec- 
trons with typical kinetic energies of the order of 1 eV. ISeatonl 
(1962) tried to rectify the Born cross-sections by modifying 
the transition probability for low-impact parameters, thus ac- 
counting for strong coupling between states, which was previ- 
ously neglected (the so-called impact parameter approximation). 
Another approach is to empirically adjust the Born rates to reach 
be tter agreement w i th exp erim ental data, w hich was advocated 
by Ivan Regemorterf (1 19621) and lPark) (119711) . Without any alter- 
natives, non-LTE applications have long had to rely almost ex- 
clusively on these simple semi-empirical formulae. 

Nowadays, much more rigorous quantum mechanical cal- 
culations can be pe rformed for simple atoms like Na. 
Igenbergs et al.1 (120081) present new convergent close-coupling 
(CCC) calculations for excitation and ionisation of neutral Na 
by electron collisions, which agree well with exp erimental data 
when available (e.g. with the measurements by iPhelps & LinL 
Il98ll) . Additionally, R-matrix calculations that are able to re- 
cover the detailed resonant structures of the cross-sections at 
threshold i mpact energies (< 5eV) have been performed by 
Gao et al ] (120101) . Here, we adopted cross-sections calculated 
with the R-matrix method in a manner similar to that described 
in Gao et al. , including seven real spectroscopic states of Nal 
(3s, 3p, 4s, 3d, 4p, 4d, and 4f) and four polarised pseudo- 
orbitals to take into account the polarisation effects between the 
scattered electron and the target electrons (see Feautrier et al. 
2010, in preparation). The calculations were performed for im- 
pact energies in the range 0-14 eV. Figure |2]compares the cross- 
sections obtained for some sample transitions with the R-matrix 



metho d and the analytical fitting functions by Igenbergs et al.l 
d2008h . The two methods agree very well; the typical difference 
between calculated rate coefficients in the temperature regime 
2000 - 20000 K is only a factor of two, although individual rates 
could differ by up to a factor of ten. 

For more highly excited states of Na I, quantum mechanical 
data are lacking and we are still dependent on general formulae. 
However, the empirically adjusted rates bv iParkJ (1197 lh perform 
surprisingly well compared with the newer data, with half of the 
transitions at 6000 K agreeing within a factor of three and 90% 
of the transitions agreeing to within a factor of ten. We adopted 
these for transitions involving levels above 5s. 

Cross-sections for ionisation by electron impact have 
also been calculated with the CCC method up to 5 s by 
llgenbergs et ail (120081) . and compare well with existing experi- 
mental data for the ground and first excited states. For ionisation 
from mor e highly ex cited levels, we relied on the general recipe 
given in lAllenl (2000). Also, to bridge the small energy gap be- 
tween the ionisation edge and the minimum energy for which 
the analytical fits to the ionisation cross-sections are valid (see 
Table 2 in Igenbergs et al. ), we assumed that th e cross-sectio ns 
for all levels were given by the general recipe dAllenl 12000b at 
the edge. 

In Sect. 13.21 we evaluate the influence from uncertainties in 
collisional rates on the statistical equilibrium of Na I by multi- 
plying and dividing all electron collisional rates by a factor of 
two, which given the present data may be indicative of the true 
uncertainty of the data for the most important transitions. 

2.4. Collisions with hydrogen 

Neutral hydrogen atoms are the most abundant atomic species in 
late-type stellar atmospheres. Still, collisional cross-sections be- 
tween hydrogen and other atoms are notorious sources of uncer- 
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Fig. 2. Three sample excitation cross-sections for collisions be- 
tween Nal and electrons. The black solid lines represent R- 
matrix calculations (Feautrier et al. 2010, Gao et al. 2010) and 
the red dashed lines represent the analytical fitting functions de- 
rived by Igenbergs et al. (2008). The transitions are marked with 
labels next to the energy thresholds. 



tainty for non-LTE appli cations, and so me empirical calibration 
of the classical recipe by iDrawinl (1 1 9681) is usually the only op- 
tion, using a general rate scaling factor Sh- Again, the situation 
for Na today is improved in this respect. 

iBelvaev et al.l (1201 Ol) have performed quantum-scattering 
calculations of cross-sections for inelastic collisions between 
H + Na with impact energie s ranging from thresh old up to 10 eV. 
Using these cross-sections, Barkle m et alJ d2010) present for the 
first time quantum mechanical rate coefficients for a large num- 
ber of bound-bound transitions in Na I, caused by collisions with 
neutral hydrogen, as well as charge exchange reactions between 
sodium and hydrog en (Na* +H Na + + H~). We use the rate co- 
efficients given in Barklem et al .1 (120101) involving all levels and 
all possible transitions below the ionic limit, i.e. up to 5p (for lev- 
els above the ionic limit, charge exchange is not possibile). The 
bound-bound rates for allowed transitions are smaller than what 
is found with the commonly used classical recipe by Drawin by 
one to six orders of magnitude, depending on the transition (see 
Fig-E). The Drawin recipe is not applicable to forbidden bound- 
bound transitions, ionisation or charge exchange. 

We estimate the bound-bound collisional rates for tran- 
sitions between hig hly excited level s with the free electron 
model described by Kaulak vsl (1 19911). which is applicable to 
Rydberg atoms (Eq. 18 in iKaulakvsl using non-hydrogenic 
wave-f unctions in momentum space ca l culate d using the meth- 
ods of iHoang Binh & Van Regemorterl dl997l)V Figure [3] com- 
pares the rate coefficients found by iBarklem et al.l yOlO) with 
those given by the Drawin recipe for transitions between low ex- 
cited states, and also compares the rates by Kaulakys and Drawin 
for transitions between highly excited states. It is easily appreci- 
ated from these figures that for Na no single scaling factor can 
be applied to the Drawin recipe to reproduce either the quantum 
mechanical rates or the Kaulakys rates. 

Barklem e t al.l (1201 Oh also provide estimates of the rate un- 
certainty, so called fluctuation factors. In Sect. l3.2l we discuss the 
effect on the non-LTE calculations when multiplying all quan- 
tum mechanical rates by the maximum fluctuation factors or by 
the suggested minimum value (0.5). 
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Fig. 3. Comparison between rate coefficients at T = 6000 K for 
collisional excitation of Nal by neutral hydrogen atoms. Only 
optically allowed transitions are shown and x-axis represents the 
energy of the transition. Top: The ratio betw een the unsealed 
Drawin formula and the free electron model of Kaul akvsl (Il99lh 
for transitions between highly excited states. Bottom: The ratio 
between the unseal ed Drawin formula an d the quantum mechan- 
ical calculations by Barklem et al.l d2010l) for transitions between 
levels up to 5s. 



2.5. Line broadening 

For line broadening caused by collis ions with neutral hydro- 
gen we use the ABO theory by lAnstee & O'Maral (l995) 
and [Barklem & O'Mara d 19971) whenever applicable, and ofh- 
er wise adopt the C(, constant for van der Waals-broadening 
bv lUnsoldl (fl955h . The latter are enhanced by a factor of 2 
(A log C(, ~ 0.8), which is typically required to match the ob - 
served line profiles in the Sun (e.g. Mas honkina et al.l |2000). 
St ark broadening constants (C4) are estima ted from the tables 
of lDimitriievic & Sahal-Brechotl(ll985lll990l) . The line data are 
summarised in Table Q]for the 11 principally considered transi- 
tions. 



3. Discussion 

3.1. Departures from LIE 

There is a great regularity in the statistical equilibrium of Na 
over the whole stellar grid covered by our analysis. As we will 
describe below, the non-LTE abundance corrections are depen- 
dent on the line strength and much of the model dependence 
can therefore be easily understood through the variation of the 
excitation and ionisation equilibrium with temperature and den- 
sity. Figure [4] shows contour plots of how the inferred cor- 
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Fig. 4. Contour diagrams illustrating how the abundance correc- 
tions vary with effective temperature and surface gravity for four 
commonly analysed lines of Na I. Only the results for £ = 2 km/s 
and [Na/Fe] = are shown. 

rections vary with effective temperature and surface gravity at 
[Na/Fe] =0. 

To describe the non-LTE line formation of Na lines, we de- 
fine the departure coefficient B x - N x /N*, which represents the 
population of a certain level x in non-LTE, divided by the corre- 
sponding population in LTE. Given a specific background con- 
tinuum opacity, the strength and shape of a spectral line are de- 
termined by the line opacity, which is proportional to the popu- 
lation of the lower level of the transition, and the source function 
S v , which is proportional to the population ratio between the up- 
per and lower leve0. Lifting the assumption of LTE, both line 
formation depth and source function may change and thus alter 
the strength and shape of the spectral line. 

The strong Nal D doublet lines at 588.9/589.5 nm are com- 
monly used in abundance studies of metal-poor stars. Indeed, in 
certain cases (warm, metal-poor dwarfs) they are the only lines 
that are sufficiently strong. The doublet originates from the res- 
onance line transition between the ground state (3s) and the two 
fine structure components of the first excited state (3pi/2 and 
3p3/2). We found empirically that as a general rule the line source 
function perfectly resembles that of a pure scattering line in a 
two-level atom, i.e. to a very good approximation Si = 1$ at 
all depths, where is the profile-averaged mean intensity. This 
merely reflects that for these lines, the photon absorption and 
emission rates strongly dominate the collisional rates and all in- 
teractions with other levels, including the continuum. Although 

5 Neglecting stimulated emission, the line source function 5; is di- 
rectly given by the expression 5//5 v = /3 I( //3; 



Fig. 5. Non-LTE abundance corrections as functions of equiv- 
alent width for selected Nal lines for a metal-poor turn-off 
star (top: T eff = 6400 K, logg = 4.0, [Fe/H] = -2, and 
£ t = 2.0 km/s) and metal-poor giant (bottom: T e ff = 5000 K, 
logg = 2.0, [Fe/H] = -2, and £ = 2.0 km/s). The black bullets 
indicate the predicted line strength and abundance corrections at 
[Na/Fe] = 



in principle for a two-level atom, Si = (1 - e Y )J$ + e v B v , the 
probability for true absorption, e v , is very close to zero at shal- 
low depths where J$ departs from B Y , so the approximation of a 
pure scattering line still holds at all depths. 

The ratio between the population of the upper and lower level 
is therefore always set by the mean intensity, and, for a specific 
line strength, this is correctly established with a simple two-level 
atom. However, the actual population of the ground level, which 
governs the line opacity and typical formation depth of the line, 
is underestimated when more highly excited levels of Na I are 
neglected. As described bv lBruls et al.l (119921) . a number of high- 
excitation levels and a ladder of high-probability transitions con- 
necting these with lower excitation levels must be established 
to obtain the correct populations. This is needed because the 
first and second excited states have photoionisation thresholds 
in the ultra-violet, where the radiation field exceeds the Planck 
function (J v > B v ), and pushes the ionisation balance to overi- 
onisation of Na I. The situation is reversed for more highly ex- 
cited levels since their photoionisation thresholds lie in the near- 
infrared regime, where the radiation field is instead subthermal. 
This photoioni sation/recombinati on picture has been describe d 
previously bv lTakeda et all (120031) and lMashonkina et a l. (2000). 

We can now qualitatively understand the non-LTE formation, 
starting with the Na I resonance lines. When either line is weak, 
it will obviously have little influence on its own radiation field. 
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Therefore, J v -B v > at shallow optical depths, as is the case for 
neighbouring continuum regions. This is governed by the tem- 
perature gradient, which determines how rapidly B v decreases 
with optical depth. At these wavelengths, the gradient is steep 
enough to produce a J v - B v excess and, consequently, a line 
source function that is stronger than in LTE (J3„/Bi > 1). This 
would tend to weaken the resonance lines compared with LTE, 
but the effect is small because the lines are formed deep in the 
atmosphere, where the radiation field is close to thermal. On the 
other hand, the photoionisation balance of Na I is always shifted 
to overrecombination in the line-forming regions. The ground 
state of Nal is thus overpopulated compared with LTE (J3i > 1), 
and the combined effect is generally a moderate line strengthen- 
ing. The resulting abundance corrections for the resonance lines 
are always minor, |AA| < 0.1 dex at line strengths below 5 pm. 
However, these lines are only that weak in extremely metal-poor 
stars. 

As the line strength increases towards saturation in any given 
model, J v drops, and the line source function becomes weaker 
than in LTE. This is naturally accompanied by a decrease in 
the excitation rate. Also the ionisation rates drop with increas- 
ing abundance as the photoionising radiation field weakens. This 
enhances photorecombination from the Na II reservoir, and all 
levels of Na I are increasingly overpopulated. The combination 
of higher line opacity and a weaker source function produce sig- 
nificant line strengthening in non-LTE, with negative abundance 
corrections as a result. This behaviour is illustrated in Fig.[5]for 
a metal-poor dwarf and giant. We note that even if the statisti- 
cal equilibrium is not necessarily pushed further from LTE, the 
abundance corrections become larger and larger as the line sat- 
urates, simply because the abundance sensitivity to equivalent 
width is small in this regime. 

With further strengthening, the line enters onto the damp- 
ing part of the curve-of-growth with the development of broad 
wings. Because photons from a wider frequency range are then 
able to excite the atoms, the excitation rate actually increases 
again, lessening the overpopulation in deep layers. The abun- 
dance sensitivity to line strength also starts to increase again, 
consequently the abundance corrections, as we define them here, 
reach a minimum value when the line is fully saturated (see 
Fig.|5j and then become less negative with higher abundance, 
although these strong lines are hardly suitable for accurate abun- 
dance analysis through equivalent width measurements. Line 
profile analysis especially of the wings is more appropriate, but 
that is not addressed here. 

Even if the two-level approximation holds true only for the 
source function of the resonance lines, the subordinate transi- 
tions (with 3p as lower level) show a very similar behaviour. 
At low line strength each line has a 'plateau' of close-to- 
constant, small abundance corrections, which becomes increas- 
ingly negative as the line saturates around 20 pm. This general 
be haviour with l i ne stre ngth has been discovered previously e.g. 
by iTakeda et al.l (120031) . The strongest subordinate transitions, 
3p-3d at 818. 3/819. 4nm, follow an almost identical behaviour 
as the resonance lines, but are offset to more negative correc- 
tions. The latter is because of the shallower temperature gradient 
with continuum optical depth in the near-infrared wavelength, 
which produces a J v - B v deficiency in the line, and a source 
function that is weaker than in LTE also when the lines are very 
weak. 

Naturally, the abundance corrections at a given line strength 
are still somewhat model-dependent. For saturated lines, the cor- 
rections are more negative for hotter models and models with 
lower surface gravity, whereas the metallicity dependence seems 
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almost negligible. Fully unsaturated lines (below 5 pm) almost 
always have corrections in the range —0.1... - 0.2 dex. 

We note the very close resemblance in the non-LTE line for- 
mation of sodium and the lighter alkali atom lithi um, whose de- 
partur es fr om LTE have bee n described e.g by ICarlsson et al.l 
(119941) and|Und_etaD(|2009). There are many striking similar- 
ities between the two elements, especially in the shape of the 
abundance correction curves. Differences mainly arise from the 
higher degree of overionisation of Li, which in turn is a direct 
result of the larger photoionisation cross-section from the first 
excited state of Li I (2p), compared to Na I (3p). The abundance 
corrections at low line strengths thus tend to be somewhat higher, 
even positive, for Li (overionisation causes underpopulation of 
Li I, thus weakening the spectral lines). 

3.2. The influence of collisions 

As described in Sects. 1231 and 12.41 the statistical equilibrium of 
Na I was calculated by accounting for collisional excitation and 
ionisation by electrons and hydrogen atoms. We now discuss the 
impact on the derived abundance corrections when varying the 
strength of collisional rates. 

When multiplying/dividing all rates for collisional excitation 
and ionisation by electrons by a factor of two, the solar equiva- 
lent widths of Nal lines change systematically by typically 1%. 
This propagates to less than ~ 0.01 dex in terms of non-LTE 
abundance corrections, and is therefore negligible. A somewhat 
larger impact is seen for hotter, higher surface gravity models, 
where electrons are more abundant. Still, the non-LTE equiva- 
lent widths calculated for a T e ff = 8000 K, log(g) = 5.0, so- 
lar metallicity model are affected by only 2^1%, corresponding 
to approximately 0.02 dex for relevant lines. Giants and cooler 
dwarfs are less sensitive, because they have lower densities of 
the colliding free electrons. The Na non-LTE calculations thus 
seem robust with respect to input atomic data for collisions with 
electrons. 

Th e quantum mechanical calculations by Bar klem et al.l 
(2010) result in rate coefficients for collisional excitation by hy- 
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drogen atoms that are lower than the commonly used classi- 
cal Drawin recipe by one to six orders of magnitude. The new 
rates thus have very little influence on the statistical equilibrium 
of Na I, in the stellar parameter range that we consider here. 
However, again analogously to Li, charge exchange reactions 
turn out to be more influential. Figure |6]illustrates how the depar- 
ture coefficients of low excited levels change for a solar model, 
when including and neglecting collisions with neutral hydrogen. 

It is conceptually correct to say that higher collisional rates 
have a thermalising effect, i.e. they tend to drive the level popu- 
lation towards LTE. However, as seen for the Sun in Fig. [6] this 
is not generally true for all levels. The populations of the ground 
state and first excited states are instead pushed farther away from 
LTE in optically thin atmospheric layers, when charge exchange 
is included. This seemingly contradictory behaviour can be un- 
derstood by realising that the collisional cross-sections of 3s 
and 3p are very small, and the changes in their level population 
are rather a secondary effect, stemming from the thermalisation 
of higher excited levels. Especially 4s, whose departure coeffi- 
cients are also displayed in Fig. [6] has a large cross-section for 
charge exchange and the level becomes almost thermalised with 
the inclusion of this process. In deeper layers, where /?4 S > 1, 
the overpopulation of the level decreases, which also propagates 
to lower excited levels. In shallow layers, where (3^ < 1, the 
added neutralisation channel through charge exchange reactions 
lessens the underpopulation of 4s, but further increases the over- 
population of lower levels. In practice, higher excited levels than 
4s influence the outcome as well, but to a lesser extent. 

Even if the statistical equilibrium is indeed influenced by hy- 
drogen collisions, this need not necessarily be the case for the the 
emergent spectral lines. This is because the source functions re- 
main unchanged, and as discussed above, the effect on the level 
populations of 3s and 3p is opposite in different regimes of the 
atmosphere, so that the net effect on the line strength is small. 
The situation is different from that of Li, for which charge ex- 
change reactions always lead to more neutralisation throughout 
t he atmosphere, with significant line s trengthening as a result 
(iBarklem et al.Ll200llLind et all 120091) . 

Charge exchange is most important in cool, high surface 
gravity models, for which abundance corrections are affected 
by up to ~0.2dex, whereas the effect in solar-type stars is only 
~0.01 dex. This is a direct reflection of the ionisation equilib- 
rium of hydrogen, because H is abundant in cool, dense at- 
mospheres. We therefore tested how the uncertainty in colli- 
sional rates affects the abundance corrections for the extreme 
case of a T e g — 4000 K, \ogg = 5.0 model. When multiply- 
ing and dividing all rates with their maximum fluctuation fac- 
tors (see Sect. l2.4K the abundance corrections vary with typi- 
cally 0.04 dex for such a model. For the Sun, the corresponding 
variation is < 0.01 dex, i.e. essentially irrelevant for practical 
purposes. 

3.3. Consequences for stellar abundance analysis 

Obviously, it is always preferable to base the abundance anal- 
ysis on unsaturated spectral lines, not least because non-LTE 
abundance corrections for fully saturated lines may reach -1 dex 
in extreme cases. This means that equivalent width measure- 
ments of the resonance lines should be used only as a last 
resort for metal-poor stars, when the 818. 3/819. 4nm doublet 
lines are not detectable or too severely blended with telluric 
lines. At metallicities in the range [Fe/H] = -2.0... - 1.0, 
the 568.2/568.8 nm doublet lines are usually to be preferred, 
whereas the 615.4/616.0 are good indicators for solar metallicity 
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stars. Given the estimated uncertainties in collisional data, the 
non-LTE modelling gives abundances for weak lines that nor- 
mally can be trusted to high precision. 

Figure [7] compares our inferred corrections with those of 
Andrievs kv et all (B007I) IShi et al l (l200l . iTakeda et all d2003h 
and Mashonkina et al. (2000) for selected models within the 
grid. Overall, the agreement between the studies is never 
worse than 0.2 dex and they usually agree to within 0.1 dex. 
Considering that all studies have (to various degree) differences 
in input atomic data, atmospheric models, and the numerical 
methods, this is not surprising. Still, all statistical equilibrium 
calculations in ID are conceptually similar and a better consen- 
sus between independent studies is clearly desired to ensure that 
accurate Na abundances can be derived. In particular, we note 
the sometimes highly significant discrepancies with the non-LTE 
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calculations published bv lGratton et al.l (1 19991) . Their abundance 
corrections show an increasing trend toward positive values (up 
to +0.5 dex, depending on the line) for very low surface gravity 
stars, which is not seen in our calculations. 

While our model atom is superior to previous studies in 
that it contains more reliable atomic data for collisions, we do 
not expect that this is the primary reason behind some of the 
larger offsets compared with previous studies. As discussed in 
Sect. 13.21 the sensitivity to the detailed collisional rates is fairly 
low for the specific case of Na in the stellar parameter range 
covered and we suspect that the differen ces in the m o del at - 
mospheric structu r e have a lar ger impact. | Takeda et al.l (1200 3b . 
iMashonkina et all d2000h and 1 Andrievsky et al.l (120071) all use 
Kurucz ATLAS9 mode ls, while IShi et al.l d2004l) use the grid by 
iFuhrmann et al.l dl997l) . The atmospheric temperature gradient 
is inevitably different for models of identical stellar parameters, 
for example because of differences in the treatment of convec- 
tion, equation of state, and opacity. This is illustrated in the top 
panel of Fig.[8j which shows the difference in temperature strat- 
ification for different versions of ATLAS9 models with respect 
to a MARCS model for a metal-poor gianfl As can be seen in 
this figure, the stratification is critically dependent on the inclu- 
sion of convective overshooting in the 1 993 grid by Kurucz. The 
most up-to-dat e versions of ATLAS9 ( Castelli & Kurucz! 12004) 
and MARCS dGustafsson et all 120081) both neglect overshoot, 
and the models shown have similar, but not identical parameters 
for the mixing length, but nevertheless show differences on the 
1 00 K level. 

The middle and lower panels of Fig. [8] illustrate the effect 
on derived LTE and non-LTE Na abundances for the 819.4nm 
line, when using different atmospheric models. In both cases the 
line strengths appear, as expected, to be connected to the tem- 
peratures around log(rgoo) = 0. The higher temperatures of the 
ATLAS9 models, especially the one with overshoot, weakens 
the line compared to MARCS for a given Na abundance (i.e. a 
given line strength generally corresponds to a higher LTE abun- 
dance). However, in non-LTE the difference is lessened because 
of a higher degree of overpopulation of neutral Na in the temper- 
ature 'jump' in the overshooting model. This in turn is caused 
by the flatter temperature gradient and consequently the less ef- 
ficient overionising ultra-violet flux in this region. This compari- 
son implies that differences between ID atmospheric models can 
produce differences of the order of ~ 0.1 dex in Na abundance 
for such a metal-poor giant. 



4. Concluding remarks 

Using our calculated abundance corrections and/or non-LTE 
curves-of-growth for 1 1 important neutral Na I lines, Na abun- 
dances that are superiour to those inferred with the LTE as- 
sumption can easily be obtained for late-type dwarf and giant 
stars. The results show a low sensitivity to uncertainties in in- 
put atomic data for collisional cross-sections, but are sensitive to 
the detailed structure of the atmosphere. To minimize the influ- 
ence from possible systematic errors in the model, unsaturated 
lines are definitely to be preferred as abundance indicators. As 
always, a good test of the soundness of the modelling procedure 
is to compare Na abundances from lines of different strengths. 



6 The ATLAS9 models were retrieved 
http://kurucz.harvard.edu/grids/ with convective 



from 
over- 



shoot (gridm20), without overshoot (gridm20NOVER), and using more 
up-to-date opacities and abundances (gridm20a2ODFNEW). 



Our results satisfactorily agree with earlier non-LTE studies 
of Na I that treated the strength of hydrogen collisions as a free 
parameter, even though the new quantum mechanical rates for 
hydrogen impact excitation are considerably lower than the tra- 
ditational estimates. This low sensitivity of non-LTE abundances 
to hydrogen collisions stems from a cancellation between effects 
in different line-forming regions. 

In a forthcoming paper we will explore the impact of the 
model atmosphere on the abundance determination more closely, 
and extend our work to include 3D, hydrodynamical model at- 
mospheres, superseding the crude mixing length recipes used for 
convection in static ID models. The model atom presented here 
will then be applied to specific cases of solar-metallicity and 
metal-poor dwarfs and giants, and we will discuss the Galactic 
chemical evolution of Na. 
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Fig. 8. Top: Temperature differe nce between ATLAS9 m odels 
with respect to MARCS models dGustafsson et all l2008h . for a 
metal-poor giant (T eff = 5000 K, logg = 2.0, [Fe/H] = -2.0, 
<;t = 2.0km/s). The dotted line represents an ATLAS9 model 
with scaled solar composition and convective overshoot from 
the Kurucz 1993-grid. The dashed line has the same composi- 
tion, but without overshoot. The dashed-dotted line represents a 
0.4 dex alpha-enhanced ATLAS9 model without overshoot from 
the 2004-grid by Castelli & Kurucz. Middle: Difference be- 
tween LTE abundance as function of line strength for the same 
models with respect to MARCS. The abundances are based 
on the 819. 4nm Nal line and vary from left to right between 
[Na/Fe] = -2.0... + 2.0. Bottom: Same as the middle panel, but 
for non-LTE abundances. 



